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Abstract. We consider general, non-linear curvature perturbations on scales greater than the Hubble horizon scale by invoking 
an expansion in spatial gradients, the so-called gradient expansion. After reviewing the basic properties of the gradient 
expansion, we derive the conservation law for non-linear curvature perturbations for an isentropic fluid. We also define the 
gauge-invariant curvature perturbation under a finite shift of time-slicing, and derive the non-linear genralization of the SN 
formalism. The results obtained are straight-forward generalisations of those already proven in linear perturbation theory, and 
the equations are simple, resembling closely the first-order equations. 



INTRODUCTION 

Analysis of the WMAP first-year data [1] revealed that 
the standard inflationary universe scenario in which 
a single inflaton field slowly rolls down its potential 
hill during inflation is perfectly consistent with the ob- 
served CMB temperature spectrum [2]. Namely, the 
CMB anisotropy theoretically predicted from a scale- 
invariant adiabatic perturbation in a spatiall flat universe 
matched with the observed CMB anisotropy with impres- 
sive precision [3], and no non-Gaussianity was found in 
the statistics of the observed anisotropy [4]. This not only 
means the validity of the standard inflationary theory but 
also justifies the use of linear cosmological perturbation 
theory [5, 6, 7]. 

Nevertheless, there are some features in the observed 
CMB spectrum which might not be caused simply by 
cosmic variance but that could be a signature of non- 
standard inflation such as multi-field inflation, non-slow- 
roll inflation, braneworld inflation, etc.. Also, the level of 
the Gaussianity test is still not stringent enough that there 
could be some non-Gaussianity which are to be found in 
future experiments, say, by PLANCK [8]. To deal with 
such cases properly, it is necessary to consider nonlinear 
perturbations on superhorizon scales. 

Turning to the state of our present universe, the 
WMAP data also confirmed the existence of energy in 
the form of a cosmological constant or vacuum energy, 
now called the dark energy, and it dominates the en- 
ergy density of the present universe. Because the energy 
scale of this vacuum (dark) energy is extremely small, of 
O(10~ 3 )eV, compared to the Planck scale, 10 19 Gev, the 
confirmation of its existence was a big shock to theorists, 
especially those from the particle physics/string theory 



community. 

Recently, several authors claimed that the present 
dark energy can be explained by backreaction of super- 
Hubble scale perturbations at second order [9, 10]. 

Nonlinear dynamics of cosmological perturbations on 
superhorizon scales has been discussed already by many 
authors [11]. But in the light of these new developments 
mentioned above, it is appropriate to revist the dynamics 
of cosmological perturbations on superhorizon scales, 
and clarify what can be said and what cannot be said. 

In this report, we give a fully nonlinear formulation 
of cosmological perturbations on superhorizon scales at 
leading order in the gradient expansion, based on our 
recent paper [12]. 

GRADIENT EXPANSION 

The gradient expansion is a powerful tool to deal 
with cosmological perturbations on superhorizon scales. 
It was developed by various authors for various pur- 
poses [11]. Here we first briefly describe the essence of 
it. 

The gradient expansion assumes that spatial deriva- 
tives are always smaller than the time derivative for any 
physical quantity Q: 



In cosmological situations, we have \d t Q\ ~ HQ where 
H is the inverse of the Hubble timescale (or th e gravita- 
tional free-fall timescale) and given by H ~ y/8nGp/3. 
At each point of spacetime, we therefore have a local def- 
inition of the Hubble horizon size, given by H . Hence 



the gradient expansion is valid if the quantities are slowly 
varying in space on the Hubble horizon scale. 

At lowest order in the gradient expansion, the Ein- 
stein and matter field equations become ordinary differ- 
ential equations in time. In other words, the evolution on 
scale of each Hubble size region, which we call 'local' 
evolution, should not depend on what is happening in 
some spatially distant part of the universe, as illustrated 
in Fig. 1. This is just a consequence of causality. 



Metric 

We consider the standard (3 + 1 )-decomposition of the 
metric, which applies to any smooth spacetime: 



ds z 



-.J^ 2 dt 2 + Yijidx' + Fdt^dxi + P j dt) , (2) 



where ,jV is the lapse function, j3 ! the shift vector, and 
Yij the spatial three metric. (Greek indices will take the 
values jU , V = 0,1,2,3, Latin indices i,j = 1,2,3. The 
spatial indices are to be raised or lowered by or Yij-) 
In this (3 + l)-decomposition, the unit time-like vector 
normal to the x° = t =const. hypersurface has the 
components, 



(3) 



We write the 3-metric, Yij> as a product of two terms, 



2a - 

Yij = e Yij, 



(4) 



where a and fij are functions of the spacetime coordi- 
nates (t,x l ), and det^] = 1. 

The above form of the metric is completely general. 
Now we apply it to the case of our interest, i.e., to cos- 
mological situations in which the metric varies substan- 
tially on scales much greater than the Hubble horizon but 
is almost homogeneous and isotropic (and spatially flat) 
on scales of the Hubble horizon. Then, although it is not 
necessary, it is convenient to introduce a fiducial 'back- 
ground' scale factor a(t ) and a perturbation y/ on it as 



(5) 



Note that the variable y describes the curvature pertur- 
bation in the limit of linear theory, 



^R(Y) = -^ i3) Ay + 0(y 2 ). 



(6) 



Hence we also call it the curvature perturbation, although 
there will be no simple relation between \j/ and in the 
nonlinear case. 

Likewise, it is convenient to factorize the matrix fij as 



f = Ie H , 



(7) 



where / is the unit matrix. The condition det(f) = 1 
ensures that the matrix H is traceless. In the limit of 
linear theory, the transverse part of Hy describes the 
tensor (gravitational wave) perturbation. 

To invoke the gradient expansion, we associate an ex- 
pansion parameter e with each spatial derivative, <9, — > 
edi. Physically, e is equal to the ratio of the Hubble 
horizon size to the typical wavelength of a perturbation, 
e ~ H~ l /X. As for the form of the metric (2) is con- 
cerned, in order to apply the gradient expansion, the only 
non-trivial assumption we make is that 



j3' = 0(e). 



Matter 



(8) 



The energy-momentum tensor is assume to have the 
perfect fluid form 



T^v = (p +P)u^u v +g^ v P , 



(9) 



where p — p (x^) is the energy density and P = P(x fl ) is 
the pressure. Then the energy conservation equation is 

- MM v v r" v = ^p + (p+p)e = o, (io) 

where djdx — u^V^ and 



e = . 



(ii) 



In accordance with the assumption (8) on the shift 
vector /3 ! , we also assume that the 3-velocity of the 
matter to satisfy 



^„-0(.). 



Then it follows that 

3d t a 



= ^ r + O(e z ) = e n + O(e z ), 



where 



(12) 



(13) 



(14) 



That is, the expansion of the matter 4-velocity and that 
of the unit vector normal to the hypersurface coincide 
to each other at lowest order in the gradient expansion, 
which follows from the assumptions (8) and (12). 

One can then introduce the notion of the local Hubble 
parameter H by 



H = -e = -9„ + 0(e z )^d t a + 0(e z ). 



(15) 



3 3 

The second (or the last) equality implies that this defin- 
tion is unique in the sense that the local Hubble param- 
eter can be defined in terms of the expansion of the hy- 
persurface normal , and the definition is independent 
of the choice of the time slicing. Here, we note again that 
'local' means 'on scales of the Hubble horizon size.' 



L »H- 1 



light cone 




FIGURE 1. Schematic spacetime diagram for superhorizon scale cosmological perturbations. Each 'local' Hubble-horizon region 
evolves independently because of causality. 



Local Friedman equation 

So far, we have not specified the theory of gravity. 
Let us now invoke Einstein gravity. By explicitly writ- 
ing down the Einstein equations, one sees that one can 
further consistently assume [12] 

^- t fij = 0(e 2 ). (16) 

Physically this assumption means the absence of a de- 
caying mode in the shear of the hypersurface. In the in- 
flationary universe, this decaying mode dies out rapidly 
soon after the comoving scale of interest leaves out of the 
horizon. We then find that the only non-trivial equation 
is the (O,0)-componentof the Einstein equations, i.e., the 
Hamiltonian constraint on the t =constant hypersurface, 
which reads 

3H 2 = 87tGp + 0(e 2 ). (17) 

This is exaclty the same as the usual Friedmann equation 
for a homogeneous and isotropic, spatially flat universe. 
Together with the energy conservation (10), we now 
see that the Friedmann equation is locally valid in each 
Hubble-horizon size region even under the presence of 
nonlinear perturbations on superhorizon scales. 

An immediate consequence of Eq. (17) is that the slic- 
ing for which the energy density is uniform on each time 
slice (uniform density slicing) is equivalent to the slicing 
for which the local Hubble parameter is uniform (uni- 
form Hubble slicing), up to errors of 0(e 2 ). Although not 
apparent from Eq. (17), one can also easily show from 
the momentum constraint equations that the comoving 
slicing for which T ; = {u°) 2 (p +P)(vj + j3j) = coin- 
cides also with the uniform Hubble or uniform density 
slicing up to errors of 0(e 2 ) [12]. 

Another, most important consequence of Eq. (17) is 
that the local physics cannot be affected by the pres- 
ence of superhorizon scale perturbations, no matter how 
large they are. In particular, this implies that there will be 
no modification or backreaction to the local Friedmann 



equation due to superhorizon scale perturbations, con- 
trary to the claim made in [9, 10]. More explicit investi- 
gations of the issue have been done recently by Hwang 
and Noh [14] and Hirata and Seljak [15], whose results 
are in agreement with our conclusion. 



NONLINEAR CURVATURE 
PERTURBATION AND AN FORMLA 

Now we investigate the evolution of the curvature pertur- 
bation y/ at lowest order in the gradient expansion. We 
first show the conservation of nonlinear curvature pertur- 
bation for an isentropic (adiabatic) fluid. We then derive 
the AN formula, which relates the amplitude of curvature 
perturbation to the perturbation in the e-folding number 
of cosmic expansion, for the general case. 

Nonlinear conservation of £ 

Multiplying the energy conservation equation (10) by 
jV , we have 

a - + d lV =- 1 --^- + 0(e 2 ). (18) 
a 3 p +P 

So far, we have not specified the choice of the time- 
slicing. Let us now choose the uniform density slicing, 
p = p(t). Following the standard notation [16], we intro- 
duce 

£ = —y/ on uniform density slices. (19) 

Then, provided that P is a unique function of p (the 
'adiabatic pressure' condition), Eq. (18) shows that d t \j/ 
is spatially homogeneous to first order, 

-d,y/ = «?,C = 0(e 2 ), (20) 

where the part of y that may be a function of only time is 
assumed to be absorbed in the fiducial background scale 



factor a{t) without loss of generality. Thus £ is conserved 
on superhorizon scales. This is a generalization of con- 
served curvature perturbation on uniform density slices 
in linear [16] and second order [13] to nonlinear order. 

It may be noted that, as in linear or second order the- 
ory, the conservation law holds for each fluid for each 
corresponding uniform density slicing, if there exist mul- 
tiple fluids, as long as they are interacting only gravita- 
tionally, and it is derived solefy from the energy conser- 
vation law, independently of the gravitational theory one 
considers. 



AN formula 

Let us define the number of e-foldings of expansion 
along an integral curve of the 4-velocity (a comoving 
worldline): 



1 f'2 1 f'2 

N(t 2 ,ti;x') = - \ OjVdt = -- dt 
3 Jti 3 Jt, 



dtp 



p+P 



(21) 



where, for definiteness, we have chosen the spatial co- 
ordinates {x 1 } to be comoving with the fluid. It is very 
important to note that this definition is purely geomet- 
rical, independent of the gravitational theory one has in 
mind, and applies to any choice of time-slicing. 
From (15) we find 



y(t 2 ,x') - y{h,x l ) = N(t 2 ,h;x i )-]n 



a(t 2 ) 



a(h) 



■ (22) 



Thus we have the very general result that the change in 
y, going from one slice to another, is equal to the dif- 
ference between the actual number of e-foldings and the 
background value A^o(f2,fi) = \n[a{t 2 ) / a{t\)\. One imme- 
diate consequence of this is that the number of e-foldings 
between two time slices will be equal to the background 
value, if we choose the 'flat slicing' on which y/ = 0. 

Consider now two different time-slicings, say slicings 
A and B, which coincide at t = t \ for a given spatial point 
x' of our interest (i.e., the 3-surfaces T. A (t\) and Eg(fi) 
are tangent to each other at x'). Then the difference in the 
time-slicing at some other time t — t 2 can be described 
by the difference in the number of e-foldings. From 
Eq. (22), we have 

VA{t2,j)-VB{hj) = N A {t 2 ,h\tf)-N B {h,h\j) 
= ANab^x 1 ), (23) 

where the indices A and B denote the slices A and B, 
respectively, on which the quantities are to be evaluated. 

Now let us choose the slicing A to be such that it starts 
on a flat slice at t = t\ and ends on a uniform-density 
slice at t = t 2 , and take B to be the flat slicing all the time 



from t = t\ to t = t 2 . Then applying Eq. (23) to this case, 
we have 

\l/ A (t 2 ,x i )=N A {t 2l t 1 ;x i )-N (t 2 ,ti)=AN F (t 2 ,ti;x i ), (24) 

where ANF{t 2 ,ti;x l ) is the difference in the number of 
e-foldings (from t = t\ to t = t 2 ) between the uniform- 
density slicing and the flat slicing, as illustrated in Fig. 2. 
This is a non-linear version of the AN formula that gen- 
eralises the first-order result of Sasaki and Stewart 

Now we specialise to the case P — P(p). In this case, 
Eq. (22) reduces to 



y{t 2 ,x i )-y(h,x i ) = -\n 



a{h) 



a{h) 



1 rPfoS) dp 



1 rP('7 
3 Jp(t u 



x-) p+P 



(25) 



Thus, there is a conserved quantity, which is independent 
of the choice of time-slicing, given by 



1 fPt'S) dp 



Pit) P+P 



(26) 



In the limit of linear theory, this reduces to the conserved 
curvature perturbation in the uniform-density, uniform- 
Hubble, or the comoving slicing, 



-a x i ) = <% c (x i ) = xi/(t,x i )+ 



Sp(t,x l ) 
3 (P+P) 



(27) 



where $% c is the curvature perturbation on the comoving 
slices [17]. 



CONCLUSION 

We have investigated the behavior of nonlinear cosmo- 
logical perturbations on superhorizon scales by invoking 
the gradient expansion. Already at lowest order in the ex- 
pansion, we have obtained some non-trivial, very useful 
results. They are summarized as follows: 

• The Friedmann equation for a spatially flat universe 
holds locally (on scales of the Hubble horizon size), 
no matter how big the perturbation is on superhori- 
zon scales. A direct consequence of this is that there 
will be no backreaction effect from super-Hubble 
perturbations on local (horizon-size) physics. Local 
physics is determined solely by local physics. 

• There exists a non-linear generalization of £, which 
describes the curvature perturbation on uniform 
density slices, which is conserved for a barotropic 
fluid on super-Hubble scales. 

• There exists a non-linear generalization of A/V for- 
mula, which relates the final amplitude of the cur- 
vature perturbation to the difference in the e-folding 
number between 'flat' and uniform density slices at 
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FIGURE 2. A schematic space-time diagram for the AN formula. The slicing that interpolates between and S(?2) is the 
uniform density slicing, while the relevant slicing for the AN formula is the one that starts from the flat slice (f i ) and ends on the 
uniform density slice ^(ti)- The curvature perturbation on the final hypersurface E( f 2) is given by the difference in the e-folding 
number from ) to Efe) and that from Y.{t\ ) to Efo). 



an initial epoch (which may be chosen arbitrarily 
provided that the comoving scale of interest is be- 
yond the Hubble-horizon scale). This formula may 
be useful in evaluating non-Gaussianity from infla- 
tion. 

In this report, we only discussed the properties of non- 
linear perturbations on superhorizon scales at lowest or- 
der in the gradient expansion. However, in non-standard 
models of inflation, such as in a non-slow-roll model, 
the second order corrections in the gradient expansion 
is known to be important already in linear perturbation 
theory [18, 19]. Thus extending the present analysis to 
second order in the gradient expansion will be necessary 
to deal with more general cases. This issue is under in- 
vestigation [20]. 
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